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The bulk microwave conductivity of a dirty rf-wave superconductor is known to depend sensitively 
on the range of the disorder potential: long-range scattering enhances the conductivity, while short- 
range scattering has no effect. Here we show that the three-terminal electrical conductance of 
a normal-metal-d-wave superconductor-normal-metal junction has a dual behavior: short-range 
scattering suppresses the conductance, while long-range scattering has no effect. 

PACS numbers: 74.25.fc, 74.45. -fc, 74.62.En, 74.72.-h 



I. INTRODUCTION 



As pointed out by Lee in an influential paper [I], dis- 
order has two competing effects on the microwave con- 
ductivity of a layered superconductor with d-wave sym- 
metry of the pair potential. On the one hand, disorder 
increases the density of low-energy quasiparticle excita- 
tions, located in the Brillouin zone near the intersection 
of the Fermi surface with the nodal lines of vanishing ex- 
citation gap. On the other hand, disorder reduces the 
mobility of these nodal fermions. For short-range scat- 
tering the two effects cancel [Ij, producing a disorder in- 
dependent microwave conductivity (Tq ~ (e^//i)fcFCo per 
layer in the low-temperature, low-frequency limit (with 
^0 the coherence length and kp the Fermi wave vector). 
For long-range scattering the flrst of the two effects wins 
[3l|4j, which explains the conductivity enhancement mea- 
sured in the high-T^ cuprates 1^ (where long-range 
scattering dominates |7])- 

The microwave conductivity is a bulk property of an 
unbounded system, of length L and width W large com- 
pared to the mean free path A finite system makes it 
possible to study the crossover from diffusive to ballis- 
tic transport, as L and W become smaller than I. We 
have recently shown [5] that the transmission of nodal 
fermions over a length L in the range £,o L <^ l,W 
is pseudodiffusive: The transmission probability has the 
W/L scaling of a diffusive system, even in the absence 
of any disorder. The corresponding conductance Go is 
close the value {W/L)ao which one would expect from 
the microwave conductivity, up to a small correction of 
order (fci^fo)"^ < 1- 

It is the purpose of this paper to investigate the ef- 
fects of disorder on the pseudodiffusive conductance, as 
L becomes larger than I. We find a qualitatively different 
behavior than for the microwave conductivity, with an ex- 
ponentially suppressed conductance in the case of short- 
range scattering and an unaffected conductance G ~ Go 
for long-range scattering. 
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FIG. 1: Geometry to measure the transmission of quasipar- 
ticles at the nodes (red circles) of the pair potential with d^y 
symmetry. 



II. FORMULATION OF THE PROBLEM 

The geometry to measure the transmission of nodal 
fermions is illustrated in Fig. [T] It consists of a super- 
conducting strip S between two normal metal contacts 
Ni and N2. The transverse width W of the supercon- 
ductor is assumed to be large compare to the separation 
L of the NS interfaces, in order to avoid edge effects. 
Contact Ni is at an elevated voltage V , while S and N2 
are both grounded. The current I2 through contact N2 
measures the transmitted charge, which is carried en- 
tirely by nodal fermions if L ^ ^q. The nodal lines are 
the X and y axes, oriented at an angle a relative to the 
normal to the NS interfaces. There are four nodal points 
A, B, C, D in the Brillouin zone, at the intersection of the 
nodal lines and the Fermi surface. The nodal fermions 
have an anisotropic dispersion relation, with a velocity 
vp parallel to the nodal axis and a much smaller velocity 
wa = vp/kp^Q perpendicular to the nodal axis. 

The (three-terminal) conductance G — I2/V was cal- 
culated in Ref. [H]in the clean limit L I, with the result 
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(2.1) 



independent of a. The factors ri,r2 G (0,1) are the 
(mode-independent) transmission probabihties of tunnel 
barriers at the NiS and N2S interfaces. We have assumed 
that the tunnel barriers do not couple the nodes, which 
requires a ^> S,o/L and ■n/A — a 3> (.o/L. Since io/L <C 1, 
this is the generic case. 

We now wish to move away from the clean limit and 
include scattering by electrostatic potential fluctuations. 
We distinguish two regimes, depending on the magnitude 
of the correlation length l^. of the potential fluctuations. 
In the regime kplc S> 1 of long-range disorder, the nodes 
remain uncoupled and can be treated separately. We con- 
sider this regime of intranode scattering first, and then 
include the effects of internode scattering when l^. be- 
comes smaller than 1/kp. 



III. INTRANODE SCATTERING REGIME 

In the absence of internode scattering, the electron and 
hole components of the wave function ^ = (^fc^/i) of 
nodal fermions (at excitation energy e) are governed by 
the anisotropic Dirac equation iJ^ — eVP. Near node A 
the Hamiltonian takes the form |S] 

H = -ih{vFcrzdx + VACTxdy) + Vfj^a^ + Fact^,. (3.1) 

The two terms V^{x,y) and V\{x,y) describe, respec- 
tively, long-range disorder in the electrostatic potential 
and in the s-wave component of the pair potential. These 
two types of disorder preserve time-reversal symmetry. 
The Hamiltonian anti-commutes with the Pauli matrix 
tTj,, belonging to the chiral symmetry class AIII of Ref. 

m 

Following Refs. irUlfTTl at zero energy, the disorder po- 
tentials can be transformed out from the Dirac equation 
by means of the transformation ^ 1-^ exp(i(/) + xCTy)5'o, 
with fields (j) and x determined by 

vpdx(j> -t- v/^dyX = -V^/h, (3.2a) 
vpdxX- VAdy<j> = VA/fi- (3.2b) 

If ff* = then also i?o*o = 0, where Hq is the Dirac 
Hamiltonian without disorder (V^ = and Va = 0). 

The transformation from to leaves the particle 
current density unaffected but not the electrical current 
density: The particle current density j reads 

{jx,jy) = '^H'"F(^z,VA<^x)'^ = 'i'l{vF(^z,VA(Jx)'^0, 

(3.3) 

while for the electrical current density i one has 

iy = 0, = ewF*^* = e?;F*oexp(2xcr3,)*o- (3.4) 



This is consistent with the findings of Durst and Lee [3] , 
that the low-energy effects of intranode scattering on the 
density-of-states and on the mobility cancel for the ther- 
mal conductivity (proportional to the particle current) 
but not for the electrical conductivity (which is increased 
by disorder). 

As we now show, for the conductance of a finite sys- 
tem, the effect of intranode scattering is entirely differ- 
ent. Following Ref. O, the conductance is determined 
by the transfer matrix M relating right-moving and left- 
moving states $1 = ($^, in Ni to right-moving and 
left-moving states $2 = {^t ^ ^2) ^2- It is convenient 
to rotate the coordinate system from x and y along the 
nodal axes to coordinates s and t perpendicular and par- 
allel to the NS interfaces. The transfer matrix is defined 

by 



^2{L,t) = / dt' Mit,t')<i>i{0,t'). 



(3.5) 



For wave vectors in the normal metal coupled to node 
A, the right-movers are electrons and the left-movers 
are holes so an electron incident from contact Ni can 
only be transmitted into contact N2 as an electron, not 
as a hole. The corresponding transmission matrix t^e is 
determined by the transfer matrix via 



tee - ^yWiiJ , AA - ^^^J . 



(3.6) 



The contribution Ga to the electrical conductance from 
node A then follows from 



Ga 



2e^ 



Trf 



(3.7) 



with a factor of two to account for both spin directions. 
The full conductance contains an additional contribution 
from node B, determined by similar expressions with a 
replaced by a — 7r/2. 



The Hamiltonian (3.1 1 does not apply within a coher- 



ence length ^0 from the NS interfaces, where the deple- 
tion of the pair potential should be taken into account. 
We assume weak disorder, I ^ ^0, so that we can use 
the clean-limit results of Ref. 8 in this interface region. 
For simplicity, we do not include tunnel barriers at this 
stage (Fi = r2 = 1). The transfer matrix through the 
superconductor is then given by 

M = exp(i(?!)R -I- <7yXR) exp{-iLvFVAVa'^o-ydt + Lipadt) 
X exp(-i(?!)L - cTj^Xl), (3.8) 

with the abbreviations 



' Vp cos^ a + sm a, 
ipa = \v^'^{v], - v\)s\\i2a. 



(3.9) 
(3.10) 



The fields 0l [t] , xl {t) are evaluated at the left NS inter- 
face (s = 0) and the fields 0R(i),XR(^) are evaluated at 
the right NS interface (s = L). 
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We now follow Ref. [TU] and use the freedom to im- 
pose boundary conditions on the solution of Eq. (3.2 1. 
Demanding x = on the NS interfaces fixes both x 
and (t) (up to an additive constant). The transfer ma- 



trix (3.8) then only depends on the disorder through the 
terms exp(i(/)pj) and exp(— i(/)L), which arc unitary trans- 
formations and therefore drop out of the conductance 



(3.7). We conclude that the electrical conductance (2.1) 



is not affected by long-range disorder. 

Tunnel barriers affect the conductance in two distinct 
ways. Firstly, at both NS interfaces, we need to consider 
all four states <I>^,j that have the same component of the 
wave vector parallel to the NS interface (3>^, have the 
opposite perpendicular component than $J,<i>^). How- 
ever, only one right-moving and one left-moving super- 
position of these modes, is coupled by the transfer 
matrix to the other side of the system: 



(2 - r„)-i/2 [cf>+ + (1 - r„)i/2$+] , (3.11a) 



$,--(2-r„)-i/2[(i-r„)i/2$^ 



(3.11b) 



The superposition of incoming electron and hole states 
orthogonal to <i>+ is fully reflected by the tunnel barrier 
and the superconductor, and so plays no role in the con- 
ductance. For a detailed derivation of these formulas see 
Appendix [A] 

Secondly, the modes $^ are only partially transmitted 
through the barriers. We have calculated the transmis- 
sion probability (see Appendix [A| for details), and found 
that it can be accounted for by the following transforma- 
tion of the transfer matrix, 



M ^ e^^^^Me^^^y, 7„ = 5 ln(2/r„ - l) 



(3.12) 



With tunnel barriers, the transmission matrix contains 
mixed electron and hole elements, 



\the thhj ^ 







(3.13) 



where the unitary matrices C/„ transform from the 
electron-hole basis to the basis state and its (fully 
reflected) orthogonal complement. 



^1/2 



(1 



1 



^1/2 ^ ) • (3.14) 



Finally, the contribution Ga to the electrical conductance 
from node A follows from 



2e2 



(3.15) 



With tunnel barriers, not just nodes A and -B, but nodes 
G and D also contribute to the full conductance. 



Collecting results, we substit ute Eq . (3.8) (with xl and 
Xr both fixed at zero) into Eq. ( 3.12 ) to obtain the trans- 
fer matrix, and then substitute the 1, 1 block into Eq. 
(3.13) for the transmission matrix. Disorder only en- 
ters through the factors exp{i(j)^) and exp(—i(j)i^) , which 



mix the modes on the superconducting side of the tunnel 
barriers. Since the tunnel probabilities are assumed to 
be mode independent, these factors commute with the 
UnS and cancel upon taking the trace in Eq. (3.15). We 
thus recover the clean- limit result (2.1), independent of 
any disorder potential. Disorder would have an effect on 
the conductance for mode-dependent tunnel probabili- 
ties, but since the modes in the normal metal couple to 
a narrow range of transverse wave vectors in the super- 
conductor, the assumption of mode-independence is well 
justified. 

As an aside we mention that the thermal (rather 
than electrical) conductance Gthermai oc Tr TT^ would 
be independent of disorder also for the case of mode- 
dependent tunnel probabilities, since the C/„'s drop out 
of the trace. The tunnel barriers would then still en- 
ter in the transfer matrix through the terms e'^""'" in Eq. 
( 3.12 ), but these terms have the same effect as delta func- 



tion contributions to and can therefore be removed 



by including them in Eq. (3.2). The conclusion is that 



the thermal conductance is independent of both disor- 
der and tunnel barriers, while the electrical conductance 
is independent of disorder but dependent on tunnel bar- 
riers through the factors r„/(2 — r„). Notice that the 
Wiedemann-Franz relation between thermal and electri- 
cal conductance does not apply. 



IV. EFFECT OF INTERNODE SCATTERING 

So far we have only considered intranode scattering. 
For short-range disorder we have to include also the ef- 
fects of internode scattering. Internode scattering sup- 
presses the electrical conductance, measured between the 
normal metals Ni and N2, because an electron injected 
from Ni into nodes A or B and then scattered to nodes 
G or D will exit into N2 as a hole, of opposite electrical 
charge. (The charge deficit is drained to ground via the 
superconductor.) The thermal conductance, in contrast, 
remains unaffected by internode scattering because elec- 
trons and holes transport the same amount of energy. 
(Again, the Wiedemann-Franz relation does not apply.) 

We first give a semiclassical analytical theory, and then 
a fully quantum mechanical numerical treatment. 



A. Semiclassical theory 

We assume that the mean free path I for intranode 
scattering is short compared to the internode scattering 
length. Semiclassically we may then describe the inter- 
node scattering by a (stationary) reaction-diffusion equa- 
tion for the carrier densities n^, 

)=0. (4.1) 

The labels v' € {A, B, G, D} indicate the nodes, with 
diffusion tensor D,y and scattering rate j^^i from to v. 
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For simplicity we assume there is no tunnel barrier at the 
NS interfaces, and seek a solution n^{s) with boundary 
conditions 

nu{0)^liS,.A + S..,B)eVpF, n,{L) = 0. (4.2) 

Here pf is the density of states per node at the Fermi 
energy, and we have chosen the sign of the applied voltage 
V such that electrons (rather than holes) are injected into 
the superconductor from Ni. 

The diffusion tensor is diagonal in the x — y basis, with 
components Z?^ and in the direction ofv^ and v/^^ re- 
spectively. The average diffusion constant is D = ^ {D^ + 
Da) and we also define Da — cos^ a + sin^ a. We 
distinguish internode scattering between opposite nodes, 
with rate 71, and between adjacent nodes, with rate 72. 
Because the solution n^{s) in the s — t basis is indepen- 
dent of the transverse coordinate t, we may replace the 
Laplacian V • Di, ■ V D^cP /ds^ with Da = Dc = Da 
and Db = Dd = 2D - Da- 

We seek the current into N2, given by 

I2 = -eW lim IDaua + Dbub - Dene - Ddud] ■ 

(4-3) 

This can be obtained by integrating the reaction-diffusion 
equation (4.1 1 in the way explained in Ref.il2i The result 



/2 = eWppW 



1 



+ 



2 



V2(7i+72)D^ 
sinh V2L2(7i-f 72)/i?„ 



j2)i2D-Da) 



sinh ^2L^-fi+-f2)/{2D^Da) 



(4.4) 



In the small-L limit (when intervalley scattering can 
be neglected) we recover an a-independent conductance 
I2/V — > e'^ppDW/L, consistent with the expected result 
(2.1). For large L the conductance decays exponentially 
(X e"-'"/'"*'", with 



mm{Da,2D ~ Da)/{-n + 12) 



(4.5) 



the internode scattering length. For weak disorder 
{kpl ^ 1) this decay length is much shorter than the 
Anderson localization length ~ /e'^^', so we are justified 
in treating the transport semiclassically by a diffusion 
equation. 



B. Fully quantum mechanical solution 

The Hamiltonian in the presence of internode scatter- 
ing belongs to symmetry class CI of Ref. [S], restricted 
by time-reversal symmetry and electron-hole symmetry 
— but without the chiral symmetry that exists in the 
absence of internode scattering. 

To write the Hamiltonian H of the four coupled nodes 
in a compact form we use three sets of Pauli matrices: 



For each i = x,y, z the 2x2 Pauli matrix ai couples 
electrons and holes, % couples opposite nodes {A to C 
and B to D), and couples adjacent nodes {A to B and 
C to D). The requirements of time-reversal symmetry 
and electron-hole symmetry are given, respectively, by 

-f^n*-f^^n, {7^®<yy)n*{j.,<E)ay) = -n. (4.6) 

In the absence of disorder, the Hamiltonian is given by 

"Hcloan = Px {vfT+ (g) (Tz + VAT- (E) (7x) "S) 

+ Py{vFT-(E)az+VAT+(E)aa;)(E)jz- (4.7) 

The momentum operator is p = —ihd/dr and we have 
defined t± = ^(to ± r^), with tq the 2x2 unit matrix. 

Since the effects of disorder in the electrostatic poten- 
tial Vf^{r) and in the pair potential Va(»") are equiva- 
lent 12], we restrict ourselves to the former. The relevant 
Fourier components of V^{r) can be represented by the 
expansion 

Vf,ir) =/io(r) 



Ai3(T-)e 



- p.A{r)e 



i(kc — kB) 



(4., 



where kx is the wave vector of node X — A,B,C,D (see 
Fig.fTl). The Fourier amplitudes Ppir) are all slowly vary- 
ing imrctions of r, with correlation length ^ ^ l/kp- The 
amplitude po is responsible for intranode scattering, aris- 
ing from spatial Fourier components of V{r) with wave 
vector <^kp (long-range scattering). The other four am- 
plitudes arise from Fourier components with wave vector 
^ kp (short-range scattering). Of these internode scat- 
tering potentials, pi,p2 scatter between opposite nodes 
and P3,P4 scatter between adjacent nodes. 

The Hamiltonian H = "Hcican + "Hdisordor contains an 
electrostatic disorder contribution "Hdisordor oc cTz- Six 
combinations of Pauli matrices are allowed by the sym- 
metry (4.6), five of which have independent amplitudes: 



4 

isordcr ^ ^ ^ 
n=0 



with 



(4.9) 



Ho = Pair) [r+ «) 70 + r_ (g) 70] po{r)TQ (g) 70, 

Hi = pi{r)T+(»-f^, n2 = P2{r)T- (S)-fx, 

"^3 = l^3{r)T^ ® 7o, "^4 = Pi{r)Tx ® ix- (4.10) 

We have solved the quantum mechanical scattering 
problem of the four coupled Dirac Hamiltonians numeri- 
cally, by discretizing "H on a grid. Since the electrostatic 
potential appears in the form of a vector potential in 
the Dirac Hamiltonian, in our numerical discretization 
we are faced with a notorious problem from the theory 
of lattice fermions: How to avoid fermion doubling while 
preserving gauge invariance |13| . The transfer matrix dis- 
cretization method we use, from Ref. [M] , satisfies gauge 
invariance only in the continuum limit. We ensure that 
we have reached that limit, by reducing the mesh size of 
the grid until the results have converged. 



5 




FIG. 2: Differential conductance as a function of sample 
length, calculated numerically from the four coupled Dirac 
Hamiltonians of nodal fermions. The solid curves are at zero 
voltage and the dashed curves at nonzero voltage. If only 
intranode scattering is present (upper curves), the differen- 
tial conductance is close to the value Gciean from Eq. (2.1 1. 



Including also internode scattering (lower curves) causes the 
conductance to decay strongly below Gcioan. 



We fixed the width of the d-wave strip at W — 150^, 
oriented at an angle a = 7r/8 with the nodal lines, and in- 
creased L at fixed ^. We set the anisotropy at Vf/va = 2 
and did not include tunnel barriers for simplicity. All five 
amplitudes pLp{r) are taken as independently fluctuating 
Gaussian fields, with the same correlation length ^. The 
Gaussian fields have zero ensemble average, (/ip(r)) = 0, 
and second moment 



(hvp) ^ dr {fip{0)np{r)). 



(4.11) 



We took Ka = 1 and either Ki = K2 = K3 = = 
(only intranode scattering) or Ki = K2 — — K4 — 
0.4 (both intranode and internode scattering). The re- 
sults in Fig. [2] give the differential conductance dl2/dV, 
both at zero voltage and at a voltage oi V = 0.2 hvp/eS,. 

Without internode scattering, we recover precisely the 
analytical result dl2/dV = Gcican a,t V — 0. At 
nonzero voltages, dl2/dV rises above Gcican with in- 
creasing L, consistent with the expectations [10] for the 
crossover from pseudo-diffusive to ballistic conduction 
at y ~ hvp/eL. Internode scattering causes dl2/dV 
to drop strongly below Gcican with increasing L, both 
at zero and at nonzero voltages. The decay is approx- 
imately exponential, consistent with our semiclassical 
theory (although the range accessible numerically is not 
large enough to accurately extract a decay rate). 



V. CONCLUSION 

In summary, we have shown that the effect of disorder 
on the electrical current transmitted through a normal- 
mctal-rf-wave-superconductor-normal metal junction is 
strikingly different depending on the range of the disor- 
der potential: Long-range scattering has no effect, while 
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FIG. 3: Sketch of the normal-superconducting interface, with 
the plane wave modes taking part in conduction with a fixed 
energy and transverse momentum. To define the modes (^^^~ , 
a piece of normal metal with length — > is inserted between 
the tunnel barrier I and the superconductor S. 



short-range scattering suppresses the current exponen- 
tially. This behavior is dual to what is known [21 E] for 
the bulk conductivity, which is unaffected by short-range 
scattering and increased by long-range scattering. Be- 
cause of the exponential sensitivity c>c e^^/''"'", we pro- 
pose the setup of Fig.[l]as a way to measure the internode 
scattering length /inter- 

As a direction for future research, it would be inter- 
esting to study the transmission in the geometry of Fig. 
[1] of low-energy excitations that are not located near the 
nodal points of the pair potential. A mechanism for the 
formation of non-nodal zero-energy states in d-wave su- 
perconductors has been studied in Refs. [131 ES] • 
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Appendix A: Tunnel barrier at the NS interface 

We consider a tunnel barrier between the normal 
metal contact Ni and the superconductor. To be spe- 
cific, we describe the left end of our setup, the deriva- 
tions for the right contact follow analogously. We in- 
troduce an additional normal metal of zero length be- 
tween the tunnel barrier and the superconductor, as il- 
lustrated in Fig. [Sj For simplicity we assume translation 
invariance along the NS interface holds: then the en- 
ergy and the wave number along the NS interface are 
good quantum numbers. The tunnel barrier mixes the 4 
modes with these constants in the normal lead Ni: 
for right-/left-propagating electrons, and for right- 
/left-propagating holes, with the 4 modes with these con- 
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slants in N'j^: 



and 



We have 
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\<t>Hj 



(Al) 



Here t = yTTe''^ and t' = \/Tie'^^ are the electron trans- 
mission amphtudes, x, x' G I^, and r and r' are the elec- 
tron reflection amplitudes. 

Since the angle a between the normal to the NS inter- 
face and the nodal line is taken to be generic, ^ a <C 
— 7r/4, the modes 0^ and (f)~ cannot propagate in the su- 
perconductor. They are localized near the NS interface, 
and follow Andreev reflection: (f)~ = —i4>^- Using this, 
we can write the scattering matrix S representing the 
combined effect of the tunnel barrier and the Andreev 
reflections on the propagating modes as 



'x -t'tt' -i'tv'*' 
\ r' t'* 
t -x'iV -r'ir'* 



(A2) 



Now there are two incoming propagating modes from 
the left, but only one outgoing propagating mode to 
the right. This implies that there is a superposition 
of $+ and $^ that is reflected with unit probability 
into a superposition of $^ and . Orthogonal to 



these uncoupled superpositions are the relevant modes 

= Ue^t + '>^h^t ^iid ^1 = "^e^e + ^h^h' which are 
coupled to the propagating mode s in the superconductor. 
We can find them from Eq. ( A2 ) by just observing what 
5t and S take (0, 0, l)t to: 



1 



Af W e 



I* piX j ' 



1 / ir'*e*-^ 



, (A3) 



where Af = \/2 — Ti is a normalizing factor. For our 
setup, all phase factors here can be absorbed into the 
definitions of the plane wave modes in contact Ni, and 
we obtain Eqs. (3.11). 



Acting with S on (u*,u^,0)^ allows us to infer the 
transmission and reflection amplitudes of the relevant 
modes, from which we can obtain the transfer matrix. 



- ; I 



Ml = 



v/ri(2-ri) 



This transfer matrix can be written in a succint form with 
a real parameter 71 characterizing the tunnel barrier: 



Ml = exp[7if7y]; 71 = ^ 



1, 2-ri 



(A5) 



This and the analagous calculation for the right edge of 
the system lead directly to Eq. (3.12). 
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